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Abstract 

We investigate the time evolution of the entropy for a paradigmatic 
conservative dynamical system, the standard map, for different values of 
its controlling parameter a. When the phase space is sufficiently "chaotic" 
(i.e., for large \a\), we reproduce previous results. For small values of \a\, 
when the phase space becomes an intricate structure with the coexistence 
of chaotic and regular regions, an anomalous regime emerges. We charac- 
terize this anomalous regime with the generalized nonextensive entropy, 
and we observe that for values of a approaching zero, it lasts for an in- 
creasingly large time. This scenario displays a striking analogy with recent 
observations made in isolated classical long-range iV-body Hamiltonians, 
where, for a large class of initial conditions, a metastable state (whose 
duration diverges with f /N — > 0) is observed before it crosses over to the 
usual, Boltzmann-Gibbs regime. 

PACS numbers: 05.20.-y, 05.45.-a, 05.70.Ce 

In his critical remarks about the domain of validity of Boltzmann principle, 
Einstein stressed JlJ that the basis of statistical mechanics lies on dynamics. 
During the past century an impressive amount of work has been made to clarify 
this connection, but certainly this crucial point is far from being well under- 
stood. In this paper, following recent observations on dissipative maps, we ana- 
lyze numerically the entropy production of a far-from-equilibrium conservative 
dynamical system, and we exhibit a new scenario where anomalous effects arise, 
and that may serve as a basis for extending the usual statistical-mechanical 
formalism. 

In low-dimensional dynamical systems, both dissipative and conservative, 
intensive effort has been done in analyzing the connection between the prop- 
erty of mixing of the system and the evolution of the statistical entropy (see, 
e.g., ||, HI [|, ||, ||, 0, H and references therein). Particularly, an interesting 
observation was made in where it was found a simple connection between 
the Kolmogorov-Sinai entropy rate (the one that stems from the properties of 
mixing of the system), and the statistical entropy (functional of a probability 
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distribution). It was shown in fact that, partitioning the phase space in W cells 
and starting from M » 1 points within one cell (details are given later), the 
time dependence of the usual, Boltzmann-Gibbs (BG), statistical entropy 



w 

Si = -^Pilnpi, (1) 

i=l 

S (t) 

includes a linear stage whose slope K\ = lim^oc lim^i/^oo limjvf->oo } coin- 
cides with the Kolmogorov-Sinai entropy rate, calculated for instance via Pesin 
equality using the Lyapunov exponents. Nevertheless, in the context of one- 
dimensional dissipative systems (e.g., the logistic map), it was found in that 
at the edge of chaos the usual entropy (|l|) fails to exhibit with nonvanishing 
slope such a behavior (which is of course expected since the Lyapunov exponent 
vanishes). Instead, the nonextensive entropy |j| (for a recent review, see JlQ| ) 

S q = ^ i ^ eM )< ( 2 ) 

succeeds for a specific value of q < 1 (q* ~ 0.2445 for the logistic map). We 
remind that the entropic form (j^) reduces to ([!]) in the limit q — > 1. 

In Hamiltonian dynamics, one of the most studied models is the standard 
map (see, for instance, [|TT|), also referred to as the kicked-rotator model: 



x t+ i = yt + 7T- sm(2nxt) + x t (mod 1), 

Z7T 



Vt+i = Vt + 7T- sin(27rx t ) (modi), 

Z7T 



(3) 



where a € R (t = 0,1,2, ...). The system is intcgrable when a = 0, while, for 
large-enough values of a (typically \a\ > 7), it is strongly chaotic^]; in Fig. 1 we 
illustrate the phase portrait of the map (g) for different values of a. It can be 
easily verified that, for all values of a, this two-dimensional map is conservative; 
it has a pair of (cc, y)-dependent Lyapunov exponents, which only differ in sign 
(simplectic structure) and globally vanish when a vanishes. Inspired by the 
results obtained in in this letter we study, for the standard map, the time 

evolution of both the extensive (q — 1) and the nonextensive (q ^ 1) statistical 
entropies, focusing our attention on small values of the parameter a, i.e., on 
those situations where the border between the chaotic and the regular regions 
has a relevant influence. For smaller and smaller values of \a\ (say below 7), 
one encounters in fact an increasingly rich fractal-like structure, characteristic 
of a large class of Hamiltonian systems. Typically, chains of regular islands 
corresponding to elliptic points in resonance condition begin to appear, and 
between each couple of elliptic points of a chain there is an hyperbolic point, 
thus forming another chain of hyperbolic points responsible for the chaotic areas. 
Around each island there is another chain of islands of higher order and, once 
again, another chain of hyperbolic points between the islands (see, for instance, 
[ p"3[ for details). This structure is usually referred to as "islands-around- islands" . 

x We are not interested here in the accelerator-mode islands that appear at various critical 
values of a (see, for instance, Jl^ ] and references therein). 
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Figure 1: Phase portrait of the standard map for typical values of a. M — 20 x 20 
orbits (black dots) were started with a uniform distribution in the unit square and 
traced for < t < 200. 



If we think that in some sense the islands-around-islands structure plays, for a 
Hamiltonian system, the role that edge-of-chaos plays for a dissipative system, 
we verify that for the standard map its relative area of influence increases when 
a approaches 0. 

To study numerically the time evolution of the entropy, we introduce a 
coarse-graining partition of the mapping phase space by dividing it in W cells 
of equal size, and we set many copies of the system (M points) in a far-from- 
cquilibrium situation putting all the M points randomly or uniformly distributed 
inside a single cell. The occupation number Mi of each cell i Q3i=i % = M) 
provides a probability distribution Pi(t) = Mi(t)/M, hence an entropy value 
S q (t) (Vg). Using then the dynamic equations @, at each step the points 
spread in the phase space causing the entropy value to change. In order to 
obtain a consistent definition of the probabilities {pi} for all t we are analyzing, 
we impose the constraint M » W max , where W max < W is the maximum 
number of cells occupied. Finally, a global quantity over the whole (x,y) phase 
space is obtained averaging many histories. Each history is characterized by 
the choice of the cell with which we start at t = (sec |Q for a preliminary 
exhibition of this phenomenon). These cells are randomly chosen in the unit 
square. We stress here that in this way we are implementing a Gibbsian-like 
viewpoint in the sense that we obtain an ensemble description of how the sys- 
tem evolves towards equilibrium. Given a, the result of this analysis is then a 
single curve of the entropy versus time, for each entropic form S q . The initial 
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Figure 2: Schematic representation of S q (t) for different values of q, in the case of a 
strongly chaotic phase space. Only for q = 1 we have, in the limit W — > oo, a finite 
entropy production K\ (K q = for q > 1 and K q — > oo for g < 1). 

part of the curve characterizes the entropy production when the system is very 
far from equilibrium. Then comes an intermediate regime whose duration in- 
creases with W . Finally, the entropy starts approaching its equilibrium, where 
the entropy saturates at a value equal or below (W 1 ~ g — 1)/(1 — q) (value of S q 
at equiprobability) . With the aim of mathematically characterizing, for fixed 
a, the properties of the entropy production in the limit W — > oo, we perform 
the calculation for increasingly large values of W. The curves asymptotically 
approach the W — > oo curve (see Fig. 2). Numerically, the crucial point is 
to execute many times the calculation with large W satisfying the constraint 
M >> W m ax- This exhausts the memory capacity of a single, standard com- 
puter if we wish to attain W larger than say 10 6 ; also, a lot of machine-time is 
required. To overcome these problems, we implemented a distributed computing 
technique which allows for the simultaneous utilization of many computers 
linked within a network. 

Let us start with a = 7, where the phase space is characterized by a con- 
nected chaotic sea (Fig. 1(a)). In this case there is just one value of the entropic 
index (q = 1) that yields a linear time evolution of the entropy S q (t), after a 
possible transient that depends on the details of the region where the initial 
data are set, and before the effects of saturation. After the transient and be- 
fore saturation, S q (t) is a convex function for q < 1 and a concave function for 
q > 1. We schematize this analysis in Fig. 2, neglecting the initial transient. 
Moreover, as said before, during the linear stage, the slope of the BG entropy 
Si(t) coincides with the positive Lyapunov exponent. 

If we now pass to a = 5, the phase space presents two large (and infinitely 
many small) islands and a connected chaotic sea (Fig. 1(b)). Depending on 
where we set the initial data, we may now face a strongly chaotic, or a regular, 
or even an islands-around-islands region. Our averaging procedure does not 
privilege any of these regions and allows the predominant ones to emerge. The 
key point is to perform many different histories so that the average stabilizes on 
a definite curve. The resulting entropy curves reflect two different phenomena: 



4 



5 10 




5 10 r 15 5 10 *15 




5 10 M5 5 10 r 15 



Figure 3: S q (t) for a — 5 and different values of q (1000 histories were averaged). 
Full (empty) circles correspond to M/10 = W = 708 x 708 (M/10 = W = 224 x 224). 
The lines are guides to the eye. The slope of Si(t) between ti = 2 and t% — 6 is 
Ki(5)~0.97. 
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Figure 4: Si(t) for different vaiues of a < 2 (5000 to 7000 histories were averaged). Fuli 
circies correspond to M = W = 1000 x 1000 for a = 2, 1.5, 1.1 (M = W = 5000 x 5000 
for a = 0.65); empty circles correspond to M = W = 448 x 448 for a = 2, 1.5, 1.1 
(M — W — 2236 x 2236 for a = 0.65). Inset: Slopes of Si(t) (ifi(a)) and of 5 . 3 (t) 
(-K0.3 ~ 0.71 + 0.30 a s ' 6 ) in their linear regimes (see also Fig. 5). The lines are guides 
to the eye. 

the sensitivity to the initial conditions of each area and the relative extension 
of that area with respect to the whole phase space (unit square). In spite of 
the presence of large islands, the results we obtain in this case qualitatively 
coincide with Fig. 2; the only difference with the case a = 7 is that the slope 
of the linear stage of Si(t), Ki(a), is smaller, consistently with the observation 
that the positive Lyapunov exponent is approaching 0. In Fig. 3 we present the 
actual numerical results for a = 5; our estimation of K\(b) ~ 0.97 is slightly 
smaller from that found in Q because here we average over the entire phase 
space, including the islands. 

What happens if we further reduce a, approaching situations similar to Fig. 
1(c)? Let us first concentrate on the BG entropy Si(t). Fig. 4 displays that, 
for a approaching zero, the part which for large values of a was linear starts 
bending downwards. Then, after a certain characteristic time (noted t cross (aj), 
linearity is obtained once again. If we now change q, we find that only for a 
special value q* ~ 0.3, S q *(t) grows linearly for t t cross (a). In Fig. 5 we plot 
So.a(t) for a approaching zero. For small values of a and for t < t cross , the 
dynamical system appears to be trapped in islands-around-islands regions, so 
that the usual entropy, Si, is forced to grow slowly. For t above t cross , the 
normal dynamical behavior is restored; in other words, there is a crossover from 
an "anomalous" regime, where S q * grows linearly, to the "usual" regime, where, 
in accordance with the observation in S, it is Si that grows linearly (before 
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Figure 5: S . 3 (t) for: a = 1.5 (a); a = 1.1 (b); a = 0.9 (c); a = 0.7 (d); a = 0.65 
(e) (5000 to 7000 histories were averaged). Full circles correspond to M = W = 
1000 x 1000 for (a), (b), (c), M = W = 5000 x 5000 for (d), (e); empty circles 
correspond to M = W = 448 x 448 for (a), (b), (c), M = W = 2236 x 2236 for (d), (e). 
Inset: t cr oss{a) defined as the intersection of the linear part (before it starts bending) 
and a standard extrapolation of the bended part of the curves Sa.z(f)- Notice that 
t C ross{cL) diverges for a —* 0. The lines are guides to the eye. Our results suggest 
limt_>oo lim a ^o limw->oo liniM^oo ^ ~ 0.71 for q — q* ~ 0.3, whereas this limit 
vanishes (diverges) for q > q* (q < q*). 

saturation). It is relevant to notice that t cross diverges for a — > (see inset of 
Fig. 5). In order to clarify these complex phenomena, we have schematized them 
in Fig. 6. In the "usual" regime a connection with an exponential sensitivity to 
the initial condition is known; the connection of the "anomalous" regime to a 
power-law sensitivity has been preliminarily verified and is now under detailed 
investigation @ (see also §, |, |, §). 

Summarizing, we have studied the production of entropy for a well known 
low-dimensional conservative map controlled by the parameter a. It appears 
that, for fixed a ^ and in the limit t — > oo (to be in all cases taken after 
the W — > oo limit), it is the standard, BG entropy S\, which is associated with 
a finite entropy production. However, during the time interval preceding this 
extreme limit, an increasingly large interval emerges for which it is S q *, with 
q* ~ 0.3, which is associated with a finite entropy production. Consistently, 
our results strongly suggest that, in the limt^oo Hm a ^o ordering, the only 
linear regime in fact observed is that corresponding to q* ~ 0.3, whereas in the 
lim a ^o liuit^oc ordering the linear regime observed corresponds to q = 1. This 
fact provides a very appealing scenario for (meta) equilibrium thermostatistics in 
long-range Hamiltonians such as that described in |l7], [l8|]. For these systems, a 
longstanding metastable state can exist (preceding the BG one) whose duration 
diverges with N — > oo, and whose distribution of velocities is not Maxwellian, 
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Figure 6: Schematic representation of S q (t) for different values of q, in the limit 
W — ► oo, when there is a crossover between two different regimes. 

but rather a power-law. The role played by a in our present simple map and 
that played by 1/N in such long-range-interacting many-body systems, appear 
to be very similar, thus providing a dynamical basis for nonextensive statistical 
mechanics ||. 
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